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Abstract 

We consider the possibility of generation of the seeds of primordial magnetic field on 
inflation and show that the effect of the back reaction of this field can be very important. 
Assuming that back reaction does not spoil inflation we find a rather strong restriction on 
the amplitude of the primordial seeds which could be generated on inflation. Namely, this 
amplitude recalculated to the present epoch cannot exceed 10 _32 G in Mpc scales. This field 
seems to be too small to be amplified to the observable values by galactic dynamo mechanism. 



1 Introduction 

Astronomical observations show that all celestial bodies carry magnetic fields. From planets to 
interstellar medium, fields of varying strength and extension have been measured. A particular 
interesting case are galaxies, galaxy clusters and beyond, the intergalactic medium (IGM) and the 
Universe at large. These fields are of order of a few micro Gauss and they extend over kiloparsecs 
or more. Unfortunately their structure is not always simple. Besides a constant component they 
have complex structure with varying symmetry, that shows that processing has taken place since 
their appearance. 

The origin of magnetic fields is unknown and many scenarios have been proposed to explain 
them. Until recently the most accepted idea for the formation of large-scale magnetic fields was 
the exponentiation of a seed field as suggested by Zeldovich and collaborators long ago. This seed 
mechanism is known as galactic dynamo, the idea is the amplification of a tiny field created early 
enough by differential rotation of the galaxies and the subsequent generation of the galactic and 
cluster fields. 

However recent observational developments have cast serious doubts on this possibility. In fact 
there are already many reasons to believe, tough this is a possible mechanism in some cases, it 
cannot be universal pQ . Some of the reasons to think that seeding cannot be an answer are simple 
01 [3]- First, the very existence of high z galaxies with fields comparable to the Milky Way is 
incompatible with the necessary number of turns. Second, the narrowness of the distribution, most 
galaxies and clusters have fields of a few micro Gauss, and this is not compatible with the different 
number of rotations and the parameters involved in every galaxy. Furthermore, magnetic fields 
seem to increase with redshift, though the evidence is not overwhelming, the sample of Faraday 
rotations measured is now consistent with an increase and the set includes tens of galaxies showing 
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this pattern. Finally, as pointed out by Dolgov, it is difficult to create the fields in clusters since 
even the most efficient ejection from point bodies in galaxies like supernovas would have difficulty 
creating them. All put together seeding seems to be ruled out and moreover, even if the galactic 
dynamo was effective, one should justify the presence of a seed field which started the process. 

This is why the mechanism responsible for the origin of large-scale magnetic fields is looked in 
the Early Universe [3] , [I] , [5] . In this paper we will consider the generation of large-scale magnetic 
fields during inflation, which as it was noted by Turner and Widrow is a natural candidate for 
doing this job [6] . 

It is known that in the Friedmann universe the conformal vacuum is preserved if the theory 
is conformally invariant [7J. Classical elctrodynamics is conformally invariant, so that photons 
should not be produced in cosmological background. Thus the conformal invariance of the elec- 
tromagnetism must be broken to produce long wave magnetic fields via excitation of the vacuum 
fluctuations. 

Different mechanisms to break the conformal invariance of electromagnetic field were proposed 
in ref. [Hj . All of them are effectively reduced either to the appearance of the effective mass or 
time dependent coupling constant 0, [5], 0, [TO], [II], [H], [13], Q3], US' PU' OH- 

Here we will consider the broad class of models where conformal invariance is broken during 
inflation and investigate the back reaction of the generated magnetic field on the background. We 
show that this back reaction is very important and leads to rather strong bounds on the maximal 
value of the strength of primordial magnetic fields which seems not enough to explain the observed 
fields as a result of amplification of these primordial seeds by dynamo mechanism. 



2 Models 



The action for the massless vector field is 



S = -~ / F^F^^—gd'x = -] I F w F^ g o°V=gd 4 x, (1) 



4 j " » 4 

where F^ = V M /1„ — W^A^ = d^A v — d v A^ is conformally invariant. It is easy to see that 
under conformal transformation g^ v — > Sftg^v the determinant transforms as g — > Q s g and g^ v — > 
Q,~ 2 g tiv . This is the reason why in the Friedmann universe with the metric 

ds 2 = a 2 (t)) (dij 2 - 5 lk dx l dx k ) (2) 

the conformal vacuum is preserved. Therefore if we want to amplify quantum fluctuations on 
inflation and thus explain the origin of primordial magnetic fields we have to assume that either 
electromagnetic field is massive or its effective coupling is time-dependent during inflation. Both 
of these options are taken into account if we write the action in the form 

5 = J (-\l 2 F^F^ + M 2 A^ J—gd\. (3) 

Here I(t) = I(cf)(t), ...), where <fi can be the inflaton, dilaton or some other scalar field and the 
dots can be anything, for instance, invariants of the curvature (see [11], [16], [M], [H]). The 
appearance of time dependence of the coefficient in front of F 2 term is naturally interpreted as 
time-dependent coupling constant of the vector field. In fact if we write the Lagrangian density 
of the vector field coupled with a charged fermion in the standard form as 

L = -J^r + i^{d, + igA^, (4) 

where g is the coupling constant, then after rescaling the vector potential by the coupling constant 
A^ — > gA^ we bring this Lagrangian to the form 



1 = -!^^ + ^"^ + ^, (5) 
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which is "ready" for introducing a time-dependent coupling constant. Note that I is an inverse 
coupling constant and small values of I correspond to a large coupling constant g, which in turn 
would mean that we are in uncontrollable strong coupling regime. Only if / is large we can trust 
the theory. For our purposes we do not need to specify in more details the origin of the time 
dependence of / here. Note that the time-dependent effective coupling leaves the Lagrangian to 
be U (1) gauge-invariant. 

The mass term introduced by "hand" spoils gauge invariance. Only when it is generated via 
Higgs mechanism the gauge invariance is preserved. On the other hand as it was noticed already in 
[5] large enough magnetic fields can be obtained only if M 2 is negative during inflation. However, 
to generate negative mass squared term via Higgs mechanism one needs a ghost scalar field with 
negative kinetic energy [T!|], [2D]- As it is well known ghosts lead to catastrophic instabilities and 
therefore we will not exploit this possibility any further here. Instead we introduce the effective 
negative mass square terms considering the non-minimal coupling of the vector field to gravity, so 
that, 

M 2 = to 2 + £R, (6) 

where for generality we also keep "hard" mass term m 2 assuming that it is positive. 

Let us now rewrite the action (J3J) in terms of the vector potential A a — (Aq, Ai) .It is convenient 
to decompose the spatial part of the vector potential in terms of its transverse and longitudinal 
components Ai — Af + di X , where diAf — (we will be assuming summation over repeated 
indices irrespective of their position) . In the homogeneous flat universe with metric © , the action 
© then becomes 

S = \J[l 2 (ArAr + AfAAj + 2A A x '-A Q AA ~ x 'A x ') 

+ M 2 a 2 (A 2 + x A x ~AfAf)]d\ (7) 

where prime denotes derivative with respect to the conformal time rj. We will consider different 
cases separately. 

2.1 Time dependent coupling 

Let us first consider the case when M 2 = and I = I (t) . Then the variation of the action |7J 
with respect to Aq gives Aq — X ' , and the action simplifies to 

S = \ J I 2 {Af Af + AJAAJ) d± x . (8) 

Substituting the expansion 

A ^^) = E / ^(^W^^. 0) 
where ef'{k), a = 1, 2 are two orthogonal polarization vectors, into ((SJ), we obtain 

S = \T. i /M CT) (k)^ CT) (-k) (4*>UW' - fc»4*UW) *n*k. (io) 

cr=l,2 J 

Rewritten in terms of the new variable 

= ai) 

this action becomes 

s = \ E / (4 CT) vt - (k 2 ^v*>) d,d*k. (12) 
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It describes two real scalar fields with time-dependent effective masses in terms of their Fourier 
components. 

We are interested in the correlation functions of the transverse components of the vector poten- 
tial and magnetic field assuming that initially the field was in its vacuum state. The quantization 
of the fields with action (|12p is standard and we will simply summarize here the results referring 
the reader to [21 , [22 for the details. Taking into account (fTT|) and we immediately find the 
correlation function 

<0|if(r ? ,x)i-^, ? ,y)|0>=- 5 -^^/ , (13) 



where fa satisfy the equations 



(14) 

which immediately follow from action (|12| . The initial conditions for these equations corresponding 
to the initial vacuum state at r\i are 

v k } fa) = r-l—^ 4°°' fa) = n/ w fa)- (!5) 
V w fa) 

These initial conditions make sense only if lu 2 > 0. Anyway we will need them only for the short- 
wavelength modes for which lo 2 ~ k 2 . The power spectrum characterizing the typical amplitude 
squared of the invariant magnitude of the vector potential, A = \Z-AiA 1 , in the appropriate 
comoving scale A = 2n/k is 

(h s V l4 J) fal 2 fc 3 hr] 

(7=1,2 

Taking into account that the magnitude of the magnetic field is 

B 2 = -B t B l = ^F ik F ik = i [diA k diA k - d k AidiA k ) , (17) 



a 



we obtain for the power spectrum of the magnetic field 

s B (k, v) = s a (k, r,) -j = 2^ 4 ^ 2a 4 /2 ■ ( 18 ) 

cr=l,2 

that is, its amplitude decays faster by an extra power of the scale compared to the amplitude of the 
vector potential. For example, a flat spectrum for magnetic field (5b (k) = const) corresponds to 
the linearly growing towards large scales spectrum for the vector potential, that is, 5a (k, rf) cx fc _1 . 

We will need to control the back reaction of the generated electromagnetic field on the back- 
ground. With this purpose let us calculate the expectation value of the energy density equal to 
Tq component of the energy-momentum tensor: 

T° = I 2 (jF^F* - F 0a F°^j = ^ {AT Af + MM) • (19) 
Taking into account © and (fTTj) we obtain 



£EM =< |T°|0 >= -Lj J2 I 

a=l,2 



, o dk 



(20) 

Let us assume that the function / depends on time during inflation and find the resulting 
spectrum of the magnetic field at the end of inflation. For short waves with k \ r]\ 3> 1 we can neglect 
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I" /I compared to k 2 in (fT4"|l and the solution of this equation with vacuum initial conditions (1151) 
then becomes 

~-L e **<™). (21) 

Because \r)\ decreases during inflation at some moment \rjk\ — 1/fc the physical scale of the wave 
with comoving wavenumber fc begins to exceed the curvature scale and taking into account that 
k 2 <C I" /I we can write the general longwave solution of (fT4|) as 

v£\r 1 )^C 1 I + C 2 I J % (22) 

where G\ and Ci are the constants of integration which have to be fixed by matching this solution 
to (p?Tj) at 1 | ~ 1/k. Let us assume that / is a power-law function of the scale factor during 
inflation 

'-"<£)"• « 23 » 

where aj is the scale factor at the end of inflation. Taking into account that 

do 



dr] = 



Ha 2 ' 

and the Hubble constant H does not change significantly during inflation we obtain from 

^C ia n + C 2 a- n -\ (24) 

2.1.1 Strong coupling case 

In case n > —1/2 the first mode dominates and, matching solutions (f2"T|) and (|2^|) at \r]k\ ~ 1/fc, 
we find 



"""TlUJ "TkVt) ■ <25) 

where we have taken into account that at the moment rjk, when the corresponding wave crosses 
the Hubble scale, the scale factor is at — k/Hj. Substituting (l^5|) into IT5|) we obtain at the end 
of inflation 

where \ p h = af/k is the physical wavelength and Hi is the Hubble constant on inflation. This 
formula is valid for HJ 1 (af/ai) > X p n > HJ 1 , where a,i is the value of the scale factor at the 
beginning of inflation. If n = 2 the spectrum of the magnetic field is flat. For H 2 ~ 10 -12 (in 
Planck units), required by primordial inhomogeneities [21] . and If ~ O (1) , the amplitude of the 
field is the same in all scales and it is equal to 5b — 10~ 12 Planck units or ~ 10 46 G immediately 
after inflation. Later on the magnetic field is frozen and decays inversely proportional to the scale 
factor squared. To estimate how much the scale factor increases after inflation we can use the 
entropy conservation law. Assuming that inflation is followed but the dust dominated stage we 
obtain 

rrl/2 / \ 1/4 

a o 1/12-"/ ( a R 



-ff 1/u ^r - , (27) 
a f T \a f J 

where g is the number of relativistic degrees of freedom of those particles which later on transfer 
their entropy to the photons, To is the temperature of the background radiation today and an the 
scale factor at the moment of reheating. The lower bound on this ratio is obtained assuming that 
reheating happens immediately after inflation. In this case for H] ~ 1CP 6 we have oo/a/ ~ 10 29 
and correspondingly the strength of the generated magnetic field cannot exceed 10~ 12 G. 

Let us calculate the energy density of the generated magnetic field. The main contribution 
to the energy density comes from the scales exceeding HJ 1 because the contribution from the 
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subhorizon scales is renormalizcd in the leading order. In the case when the dominant mode v oc I 
and A T oc v/I oc const the time derivatives of the vector potential in (I19p contribute only in 
subleading k 2 order and their contribution is comparable to the contribution of the magnetic field 
itself given by the last term in (|2"D|) . Thus we obtain 

sem = / \v k {ri)?k A dk, (28) 

a JHja, 



where at is the value of the scale factor at the beginning of inflation. Substituting (|25p into 
we find that at the end of inflation when a = a/ 



r 5^> ™< 2 ' 

sem = O (1) Hj x J ln (%) > fOT » = 2- (29) 

/ \ 2(n-2) 

We see that the magnetic field energy can be comparable with the energy density of the back- 
ground only for n > 2. Requiring that inflation should last at least 75 e- folds we obtain that the 
contribution of the magnetic field energy density does not spoil inflation, that is, £em is smaller 
that Hj until the end of inflation, only if n — 2 < 0.2. Thus, we can have slightly growing toward 
large scales spectrum of the magnetic field. In particular, for n ~ 2.2 the amplitude of the mag- 
netic field in Mpc scales can be larger by a factor 10 5 compared to the considered above case of 
the flat spectrum, that is, 8b — 1CP 7 G today. This is the greatest amplitude of the primordial 
magnetic field which we can obtain in the considered above case. Note that the theory where I 
grows with the scale factor corresponds to the case when the effective coupling constant, which 
is inversely proportional to /, is incredibly large at the beginning of inflation and becomes of the 
order of one at the end of inflation. Hence at the beginning we are in strongly coupled regime 
where such theory is not trustable at all. 

The case considered above is the only one in which we can generate strong enough fields on 
inflation. Let us show that in all other cases there is very strong bound on the possible value of 
the generated field due to the back reaction of this field on the background. 



2.1.2 Weak coupling case 

For n < —1/2 the second term in |2"^|) dominates and 



v k {n) oc a~ n -\ (30) 

In this case the result follows immediately by substituting in the formulae ([25| and (|26|) — n — 1 
instead of n, so that 

<•» -Mir 1 - 

and 

MV ,, ;) ,_g_(_vp. (32) 

Thus the spectrum of the magnetic field is flat for n — —3. This case corresponds to the coupling 
constant growing as I^ 1 oc a 3 , that is, it changes from extremely small values at the beginning of 
inflation to values of order of unity at the end of inflation. Thus the theory is trustable everywhere. 
However here the back reaction of the field is very large because A oc v/I oc a~ 2 " -1 changes very 
fast and the main contribution to the energy density comes from the time derivative of the vector 
potential in (TT9|) . that is, from the electric field. Substituting (f3T|) in ([20]) we obtain that at the 
end of inflation 

f ^ n>-2, 

^^flfJ mfe), ^ n = -2, (33) 

87T 2 I i /n ^-2(n+2) 



-2 



n < -2. 
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Requiring that inflation should last at least 75 e-folds, we find that £em < Hj at the end of 
inflation only if n > —2.2. Thus the flat spectrum for magnetic field cannot be generated during 
inflation because in this case the back reaction of the electromagnetic field would spoil inflation too 
early. In the most favorable admissible case n ~ —2.2, the amplitude of the magnetic field decays 
as 5b oc A~°' 8 and its value cannot exceed 10~ 32 G in Mpc scales today. Thus in this model with 
weak coupling constant during inflation one cannot explain the origin of the primordial magnetic 
field. 

2.2 Massive field 

Now we set 1=1 and consider the case when magnetic fields are generated by the mass term in 
the action. Variation of action (J7]) with respect to Aq gives 

A x ' - AA + M 2 a 2 A Q = 0. (34) 

Taking the Fourier transform 

X(x,7j) = /xk(*7)e ik - x ^^, M*,V) = /^ok(?7)e tkx -^ 7 j, (35) 
we obtain from here 

A* = fc2+ fc M2fl2 X k - F kX ' k . (36) 

Substituting into action © the expansions ([H]), (|33|) and using to express ^4ok in terms x'k 
we obtain 

s = It, /(^'-(fc 2 +MV)^/>) ¥ ^ 

cr=l,2 •* 

+ 1 J Slgn (i _ Fk ) - (k 2 + M 2 a 2 - ^=J0 XkX-k^) d V d% (37) 

where is defined in (JTTJ (J = 1), and 



Xk = A^ll-FfclXk. (38) 

Thus we see that in the case of massive field the longitudinal degree of freedom \ becomes dy- 
namical. In the case of positive mass squared Fk is always smaller than unity and therefore the 
sign in front of the longitudinal part of the action is positive. However, if M 2 is negative then 
1 — Fk is negative for high momentum modes with k 2 > M 2 a 2 and these modes have negative 
kinetic energy. The low momentum modes with k 2 < M a 2 have positive kinetic energy because 
Fk is negative for them. Thus, introducing a tachyonic mass for the vector field in a "hard" 
way seems to lead inevitably to the appearance of ghost for high momentum longitudinal modes 
19J. Therefore if we want to avoid catastrophic instabilities related with ghost fields we have 
to consider tachyonic vector field only as a low energy effective field theory description of some 
unknown yet theory with "safe" ultraviolet completion. On the other hand if negative effective 
mass appears as interaction with the curvature, M 2 = £R, then the field is massless on scales 
smaller then the typical distance between particles inducing the average curvature and thus there 
is a natural ultraviolet cutoff in the theory. Note that this argument is not directly applicable 
in the presence of the cosmological constant. Let us assume that the problem of ghosts can be 
somehow solved and proceed with the calculation of the magnetic field from inflation in the theory 
with M 2 = to 2 + £R. In the case m = the photon mass is m 7 ~ i? 1 / 2 , where i? 1 / 2 ~ H. Today 
it would be m 7 = H to day ~ 10 _33 ey, well below the available experimental limits on the photon 
mass. The breaking of charge conservation also manifests itself only on scales of the horizon or 
larger (> H^ 1 ~ 10 28 cm) and hence has no observable consequences. 
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The equations of motion for transverse and longitudinal modes follow immediately from the 
action (37}: 

v£>" + (fc 2 + MV) w« = 0, (39) 

and 



Let us consider de Sitter universe where 



fc 2 + M 2 a 2 - ^EE± Uk = 0. (40) 



Taking into account that R = —12H 2 , for m 2 = equation (|39p becomes 

v?" + (fc 2 ^) = 0. (41) 

For short waves with k\rj\ 1 the solution of this equation corresponding to vacuum initial 
conditions is 

„Wfo)~^ e *<™). (42) 

For k\rf\ -C 1 we can neglect the fc 2 term in (|41| and the dominating longwavelength solution of 
this equation is 

# <■*>-&(¥)"• -5(^+ 5 «- 1 )' <«» 

where we use the matching conditions at ~ 1/k to fix the constant of integration. Since here 
the calculations are very similar to those in the previous section we can immediately write the 
result for the magnetic field 

\ \ n—2 



H 



5 B (X ph , Vf )^0(l)HU-^) . (44) 



For £ = 1/6 we have n = 1 and the spectrum linearly decays with the scale. In this case its value 
today is about 10~ 37 G in Mpc scales. The fiat spectrum is obtained for £ = 1/2. However, to find 
out whether this case is possible we have to verify that the back reaction of the magnetic field will 
not spoil inflation too early. In the energy density also contributes the longitudinal mode and to 
determine its contribution we will need a longwavelength solution for Xk- It is easy to check that 
the term which is different in the equations f3"§|) and (j4"U|) can be neglected for both shortwave and 
longwave solutions and hence 

Variation of action ([3]) , where 1=1 and M 2 = m 2 + with respect to the metric gives 

+ (m a + CiJ)A M A'' + ^^ a A o + e[*gV 0, V a (A%)-VpV(A%)]. (46) 
As a result of straightforward but rather lengthy calculations we obtain 

< 0|T °|0 >=e T + e L , 

where 

^=8^4 E / H CT) T"6£ai7|^| 2 '+(fc 2 + mV + 6^V)|^| 2 ]fc 3 ^ (47) 

u— 1,2 
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is the contribution of the transverse modes, H — a' /a 2 is the Hubble constant. The contribution 
of the longitudinal mode is given by 



e L = g^jr / (1 - F) [(1 - 6CbF) \ x ' k \ 2 - ^aH (I±| 



/ m 2 a 2 + 6fH 2 a 2 > 

l* fe l + I t^f ' lxfc 



, n dk 

(48) 



where 



X = X/VW^\, b = H + 7H * + 4H % . (49 ) 

For the longwave modes with k 2 <§C |M 2 o 2 | we have Fk <C 1, X — X an d their contribution to the 
total energy density is the same as the contribution from the transverse mode. It is interesting 
to note that the longitudinal mode is the ghost in de Sitter background. However, in Friedmann 
universe filled by matter with positive pressure it is not ghost in spite of the fact that the effective 
mass squared is negative. 

Substituting (f43|) into (|47|) we find that in the leading order the contribution of the longwave 
modes into the energy density in the case m = is 

e L ~0(l)^(n 2 -12< + 6£) f ' \v k ( V )\ 2 k 2 dk, (50) 

a JHja, 



and calculating the integral we obtain 



l 

l-n ' 



for n < 1, 



e L ~0(l)Hj(n 2 - 12^ + 6^) { 1 f a ^-D (51) 

1 ' 1 , for n > 1. 



n—X \ a 



In the case £ = 1/6 and when n = 1 the contribution is canceled in the leading order and k 2 
terms give a contribution of the order of Hj, that is the same as for n < 1. However, for £ > 1/6, 
and correspondingly n > 1, the energy density of the longwavelengh electromagnetic waves grows 
with time rather fast. It is negative and therefore when it becomes of order H 2 inflation is 
over. Requiring that inflation should last at least 75 e-fold we find that the contribution of 
electromagnetic field does not spoil inflation only if n — 1 < 0.2. Thus, in the most favorable case 
of n ~ 1.2, the amplitude of the magnetic field decays as 5b oc A~°' 8 and its value does not exceed 
10 -32 G in Mpc scales today. 



3 Conclusions 

In this paper we have studied the generation of large-scale magnetic fields in a two classes of 
models. In the first case the conformal invariance of the Maxwell field was broken by a non- 
minimal coupling of the form RA 2 , this gives a non-zero time-dependent mass to the photon. 
In the second case the conformal invariance is violated because of the time-dependent coupling 
constant, I(t)F^ v F^ Vl where I(t) — I(<fi(t), ...) is a general function of non-trivial background 
fields and <fi can be for instance inflaton or dilaton. 

In principle it looks like inflation can strongly amplify the vacuum quantum fluctuations and 
therefore can lead to sizable magnetic fields. However, if we take into account the back reaction 
of the electromagnetic field and require that inflation lasts at least 75 e-folds, the strength of the 
primordial field cannot exceed 10 _32 G on Mpc scales and it is not clear whether such a small field 
can work as a seed for a possible dynamo mechanism. 

Only in the strong coupling case, I(t)F^ u Fuv, where I = If (a/a/)™ and n ~ 2.2, the amplitude 
can reach the interesting value of 10~ 7 G today. However, this case corresponds to the situation 
when the effective coupling constant is extremely large at the beginning of inflation and becomes 
of the order of one at the end of inflation and hence the theory is not trustable. 
We conclude therefore that the models considered above are not efficient in producing primordial 
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magnetic fields during inflation and, even if the galactic dynamo was effective, the field produced 
seems to be too small to play the role of a seed for this mechanism. 
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